Formation of ultra-compact X-ray binaries through circum-binary 

disk-driven mass transfer 
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ABSTRACT 

The formation of ultra-compact X-ray binaries (UCXBs) has not been well 
understood. Previous works show that ultra-short orbital periods (< 1 hr) 
may be reached through mass transfer driven by magnetic braking in normal 
low/intermediate-mass X-ray binaries (L/IMXBs) only for an extremely small 
range of initial binary parameters, which makes it difficult to account for the 
rather large population of UCXBs. In this paper we report the calculated results 
on mass transfer processes in L/IMXBs with a circum-binary disk. We show 
that when the orbital angular momentum loss due to a circum-binary disk is in- 
cluded, ultra-short orbital periods could be reached for a relatively wide range of 
initial binary parameters. The results of our binary models suggest an alternative 
formation channel for UCXBs. 

Subject headings: binaries: close - stars:evolution - X-rays: binaries 



Introduction 



Ultra-compact X-ray binaries (UCXBs) are X-ray sources with very short orbital periods 
(P < 1 hr). They are thought to be powered by accretion from a Roche-lobe (RL) filling 
donor star to a neutron star (NS). The donor has to be compact, such as the helium (or more 
heav ier elements) core of an evolved giant star or a white dwarf (WD), to fit in the small 
RL (INelson et al.lll986l ). Spectra observations have shown possible C/0 o r He/N features 



in some of the UCXBs or UCXB candidates (iNelemans et al.l I2004J . l2006l . and references 
therein). Particular interests have been paid to UCXBs in recent y ears since they a re thought 
to be potential sources for the gravitational-wave detector LISA (jNelemand l2009l ) . 



Scenarios for the formation of the UCXBs can be summarized as follows. For UCXBs in 
globular clusters, t hey often invok e dynamical process e s, such as (1) direct col l isions between 



a NS and a giant (IVerbuntI 119871 : iDavies et al.l Il992l : iRasio &: Shapiro! Il99ll : llvanova et al. 
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2005 



a NS (IBailvn fc Grindlav 



ombardi et al.l 120061). (2 ) tidal capture of a low-mass main-sequence (MS) star by 



1987). and (3) exchange inte raction between a NS and a primor- 



dial binary (jPavies fc Hansen! Il998l : iRasio et al.l |2000| ) . All these dynamical scenarios are 
involved with the so-called common-envelop (CE) phase when the mass transfer between 
a NS and a (sub)giant is dynamically unstable, which may help form a tight NS+WD or 
NS-|-He star binary. After that, the orbital period decays to ultra-short regime through 
gravitational radiation (GR) until the WD/He star overflows its RL. For UCXBs in the 
Galactic field where dynamical collisions are not important, generally two CE phases are 
required to form a ti ght NS-I-WD or NS-l-He sta r binary, which d e cays to the ultra-short 



regime through GR (jTutukovfc YungelsonI Il993l : llben et al.l Il995l : lYungelson et al.l 12002 



Belczynski &: TaamI |2004| ). An alternative formation channel of UCXBs is through sta- 
ble mass transfer in no rmal low /intermediate- m ass X-ray binaries (L/IMXBs) driven by 
magnetic brak ing fMB: IVerbunt fc Zwaaru Il981l). which is now called the "magnetic cap- 



2002 



ture" scenario (Paczvnski fc Sienki ewiczl 1981 : Pvlvser fc Savoniie 1988 



Nelson fc Rappaportli2003. : van der Sluys et al.ll2005al ). However 



^odsiadlowski et al. 



van der Sluys et al. 



( l2005al ) found that only a very small range of initial parameters are allowed for the bi- 
naries to evolve to UCXBs through this channel within the Hubble time. This makes it 
impossible to account for the relatively large number of observed UCXBs. Another prob- 
lem of this scenario is related to the efficiency of MB. Observations of some rapid rotators 
in young, open clusters seem to be contradicted with the MB law originally suggested by 
Verbunt &: ZwaanI (119811). and a modified mag n etic braking law was propose d to resolve this 
problem ( Queloz et al.1l99"8 : Sills et al. 2000). van der Sluys et al. ( 2005b ) showed that, if 
this new magnetic braking law is adopted, no UCXBs can be formed through the magnetic 
capture channel. 

In this paper we propose an alternative scenario for the formation of UCXBs through 
mass transfer between a NS and a MS star. We assume that during Roche-lobe overfiow 
(RLOF), a small fraction of the mass lost from the donor forms a cir cum-binary (CB) disk 
around the binary rather accretes onto the NS (Ivan den Heuvellll994l ). Previous works have 
shown that a CB disk can extract orbital angular momentum from the binary effectively 
( ISpruit fc Taamll200ll : iTaam fc SpruitI 120011). and enhance the mass transfer rate, thus con- 
siderably influence the binary evolution (IChen fc Lill2006l : IChen et al.ll2006l ). In this work 
we explore the possible role of CB disks in the formation of UCXBs. In §2 we introduce 
the binary evolution code and the input physics incorporated in our calculations. In §3 we 
present the calculated results and compare them with observations. The implications of our 
model and related uncertainties are discussed in 54. 
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2. Evolutionary code and Input physics 



2.1. The steUar evolution code 



We use the STAR programme, originally developed by lEggletonl (Il97ll . Il972l ) and up- 
dated by other authors (IHan et al.lll994j : iPols et al.lll995l ). to compute the binary evolution. 
The ratio of mixing length to pressure scale height x = V-^p is set to be 2.0 and con- 
vective overshooting parameter t o be os = 0.12, imp l ying a 0.24ijp overshooting dist ance. 



The opacity table is t a ken f rom iHubbard fc Lampd (Il969l ). iRogers fc Iglesiad (Il992l ) and 
Alexander fc FergusonI (119941 ). Solar compositions {X = 0.70, Y = 0.28 and Z = 0.02) are 
adopted. The binary system is initially composed of a NS primary of mass Mi and a zero-age 
main-sequence (ZAMS) secondary of mass M2 with an orbital period Pi. The effective radius 
of the RL for the secondary star is calculated from lEggletonl (119831 ) equation. 



R 



0.49g 



-2/3, 



'L,2 



0.6g-2/3 + in(i + g-i/3)' 



where q = M2/M1 is the mass ratio of the binary components, and a is the orbital separation. 
The rate of mass transfer via RLOF is calculated with — M2 = RMT ■ max[0, {R2/Rl,2 ~ 
l)^]MQyr~^ in the code, where R2 is the radius of the secondary, and we set RMT = 10^. 



2.2. Angular momentum loss mechanisms 

We assume that the secondary star rotates synchronously with the binary orbital revo- 
lution, since the timescale of tidal synchronization is generally much shorter than the char- 
acteristic evolutionary timescale of the binaries considered here. We take into account four 
kinds of angular momentum loss mechanisms described as follows. 

(1) Gravitational radiation, which becom es important when the or bital period is short. The 



angular momentum loss rate is given by (ILandau fc Lifshitzlll975l ) 



dJ 32 MfMl{Mi + Ma)^ 



/2 



di'^"^ 5 c5 a7/2 



(2) 



where J, G and c are the orbital angular momentum, gravitational constant and speed of 
light, respectively. 

(2) CB disk. We assume that a small fraction 6{<^ 1) of the mass lost from the donor feeds 
into the CB disk rather leaves the binary, which yields a mass injection rate of the CB disk as 
Mqb = —SM2. Tidal torques are then exerted on the binary by the CB disk via gravitational 
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interaction, thus extracting the angular momentum fror n the binary system. The angular 
momentum loss rate via the CB disk is estimated to be (ISpruit &: Taanj|200ll ) 



dJ, 
dt 



CB 



-7 



27ra^ 



M, 



CB 



1/3 



(3) 



where 7^ = r\ja= 1.7 {t\ is the inner radi us of the CB disk), t is the t ime since mass transfer 
begins. In the standard a- viscosity disk (IShakura fc SunyaevI Il973l ). the viscous timescale 
tvi at the inner edge r-^ of the CB disk is given by = 2'y^P/3iiaj3'^, where a is t he viscosity 
param eter (we set a = 0.01 in the following calculations), (3 = H-Jri ~ 0.03 ( jBelle et al. 



20041 ). and Hi is the scale height of the disk. When Mcb is sufficiently large, the mass 
transfer will become dynamically unstable, so we add an ad hoc term 

exp(l + M2/2MEM) if - M2 > 2MEdd 

to Eq. ([3]) to suppress the CB disk-induced angular momentum loss rate when the mass loss 
rate is high. We find that this term does not infiuence the evolutionary tracks considerably, 
only expanding the parameter space of 5 suitable for UCXB formation within the Hubble 
time by ~ 10% — 20%. Here the Eddington accretion limit is expressed as 



MEdd = 3.6 X 10"^( 



Ml ,,0.1,, 1.7 1 



(4) 



where t] = GMi/R(? is the energy release efficiency through accretion {R the NS radius), 
and X is the mass fraction of hydrogen in the accreting material. 



(3) Mass loss. Similar as in iPodsiadlowski et al.l (120021 ). we assume that the NS accretion 
rate is limited to the Eddington accretion rate, and that when the mass transfer rate is less 
than MEdd, half of the mass is accreted by the NS, i.e.. Mi = min(MEdd, — ^2/2). The 
excess mass is lost in the vicinity of the NS through isotropic winds, carrying away the 
specific angular momentum of the NS, i.e. 



dJ ^JiMaafw, 



dt 



ML 



IM2I < 2MEdd 

\ (M2 + MEdd)a?cu, IM2I > 2MEdd 



(5) 



where oi = aM2/(Mi + M2) is the orbital radius of the NS, and u is the orbital angular 
velocity of the binary. 



(4) M agnetic braking. We use the saturated magnetic braking law suggested in ISills et al. 
f |2000l ). 



dJ, 
dt 



MB 



Jh 

Rq 



\Mq) 
\ 1/2 / \ -1/2 

Rq) \Mq) 



,UJ>UJcj: 



(6) 
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where K = 2.7 x 10^^ gcm^s ( lAndronov et al.ll2003l ). uJcr is the critical angular velocity at 
which the angular moment um loss rate reaches a saturated state, and can be estimated as 
( iKrishnamurthi et al.l 119971 ). 



cr,0" 



n 



(7) 



where cUcr,© = 2.9 x 10^ Hz, t^^^q is the global turnover timescale for the convective envelope 
of the Sun at its current age, rt for the secondary at age t, solved by integrating the invers e 
local conv ective velocity over the entir e surface convective envelope (IKim fc Demarquelll996l ). 
Following iPodsiadlowski et al.l (120021 ). an ad hoc factor is also added to Eq. ([6]) 



exp(-0.02/gcon + 1) if gcon < 0.02, 

where gcon is the mass fraction of the surface convective envelope. This term is used to reduce 
the strength of MB when the star has a very small convective envelope and hence does not 
have a strong magnetic field. 



3. Results 



3.1. Parameter space of Pj and M2,i 



Similar as in Ivan der Sluys et al.l (l2005al ). we define UCXBs as X-ray binaries with 
P < 50 min. The binary systems we considered are initially composed of a 1.4 Mq NS and a 
0.5— 5Mq ZAMS star. We have performed calculations of a large number of binary evolutions, 
to search suitable values in the three dimensional binary parameter space S (assumed to be 
constant through one evolutionary sequence). Pi and M2_i for binaries evolved to UCXBs 
within the age of the universe (13.7 Gyr). The distribution of Pi and M2,i for successful 
systems is shown in Fig. [1] with Pi ~ 0.7 — 1 d and M2_i ~ 1 — 3.5 Mq. When the donor mass 
My, \ > 3.5M0, the mass transfer beco mes dynamically uns t able (se e also IPodsiadlowski et al. 



2OO2I ) . Comparing with the results of Ivan der Sluys et al.l (l2005bl ) one can find that there is 
a relatively large parameter space for the formation of UCXBs if there is a CB disk at work. 



3.2. Limits and influence of 6 

The possible distribution of S and Pi with M2,i = I.IM0 is shown in Fig. [2l To illustrate 
the effects of 6 on the binary evolution, in Figs. Ej-El we plot the exemplarily evolutionary 
tracks of the secondary in the H-R diagram, the evolution of the donor mass and period as a 
function of age respectively, for a binary system with M2,i = I.IM©, Pi = 1.04 d and different 



- 6 - 



values of 5. The lower limit of S is determined by the constraint that age of the binary should 
be less than 13.7 Gyr: a larger value of 6 leads to shorter formation time, as seen from Fig. [S] 
if 6 is too small, the binary will not be able to reach the 50 min period within 13.7 Gyr due 
to inefficient angular momentum loss. The upper limit of 6 is determined by two conditions. 
The first is that the minimum orbital period Pmin should be less than 50 min. With larger 
values of 6, the donor will spend relatively less time on the MS, leaving a smaller degenerate 
helium core. According to the mass-radius relation of degenerate stars, the smaller mass, 
the larger radius, which corresponds to a larger Pmin- The second is 5 < 0.015, because we 
find that mass transfer becomes unstable in most cases if 5 > 0.015. 

The formation and evolutionary paths of UCXBs depend on the adopted values of 6. In 
the case of M2,i = I.IMq and Pj = 1.04 d, for example, when 6 < 0.0055, the orbital period 
ffist decreases with mass transfer until the donor star loses its outer envelope and shrinks 
rapidly at P ~ 0.1 — 0.2 d. This causes a cessation of mass transfer. In the subsequent 
evolution the orbital period may decrease down to the ultra-short regime under the effect of 
GR, until the secondary star fills its RL again, and the binary appears as a UCXB. When 
6 > 0.0055 the binary evolves directly into the ultra-short regime with decreasing orbital 
period. 

We need to mention that the distribution of S depends on the valu e of the viscous 
parameter a, which we cho s e to b e 0.01 in our calculations as adopted in ISpruit fc Taam 
( 120011 ) and iTaam fc SpruitI (I20011). This is abo ut an order of magnitude lower than the 
value (~ 0.1 — 0.4) inferred by iKing et al.l (120071 ) for fully ionized, thin accretion disks from 
observations of dw arf nova ou t bursts and outbursts of X-ray transients. However, the a 
value estimated by iKing et al.l (120071) is an average ori e over the entire disk, while here it 
is at the inner edge of t he disk (ISpruit fc TaamI 120011). which may be smaller due to the 



boundary condition (e.g. iPapaloizou fc Nelsonll2003l : [Winters et al.ll2003l : iFromang fc Nelson 



20061 ). Additionally, since the CB disk is located outside the binary, and hence shielded 
from X-ray irradiation from the NS by the accretion disk around the NS and the secondary, 
the a value in a CB disk might also be smaller than in accretion disks. Taking account 
of the above facts, we suggest to regard Fig. [2] as an optimistic case for the distribution of 
6. Nevertheless, from Eqs. (3) one can see that the CB disk-induced angular momentum 
loss rate is proportional to a^^^5. So if keeping a^^^5 constant, the binary evolution will be 
exactly the same (this has been verified by our test calculations), implying a predictable 6 
distribution for a given value of a. 
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3.3. Comparison with observations 

There are currently 10 UCXBs with known periods, 5 of which are persistent sources 
and 5 are transients. We list the orbital periods Porb and mean mass accreting rates Mi (or 
the upper limit of Mi) of these UCXBs in Table 1. The Mis listed in this table for persistent 
UCXBs are calculated by using the luminosities mentioned in the references and assuming 
accretion onto a lAMi^ N S with a 10 km radius, while those for transient sources are from the 



estimates of iKrauss et al.l (120071 ) , IWatts et al.l (120081 ) and iLasota et al.l (120081 ) . To compare 
observations with our CB disk-assisted binary model, we plot the Mi(= — M2/2) vs. Porb 
relations in Fig. [Sjfor binary systems with M2 = I.IMq, Pi = 1.04 d and 6 = 0.005 — 0.009, 
and in Fig. [7] for binary systems with M2 = I.IMq, 6 = 0.005 and Pi = 0.94 - 1.04 d. Note 
that in these two figures the values of Mi are not assumed to be limited to MEdd, so that 
possible super-Eddington accretion can be allowed. In this way we can compare our results 
with observations directly. We also indicate in Figs. [6] and [7] whether the accretion disks 
in the LMXBs are thermally and viscously stabl e, according to the sjtability criterion for a 

We use the symbols 



mixed-composition {X = 0.1, Y = 0.9) disk from lLasota et al 
X, *, and + on the evolutionary tracks to denote where the hydrogen com position X of 



the do nor becomes 0.3, 0.2, and 0.1, respectively, to show that the criterion of iLasota et al. 



(12008! ) is applicable here. The positions of UCXBs are marked in these two figures with 
circles and triangles for persistent and transient sources, respectively. B esides them, we 
also include 18 NS LMXBs with kno wn P and Mi (data are taken from iLiu et al.l l2007l : 
Watts et allbood : Ifleinke et al]l2009h . 



A comparison between our CB disk-assisted binary models and the observations of (com- 
pact) NS LMXBs suggesting that it is possible to form UCXBs from normal LMXBs. We note 
that three of the UCXBs are in globular clusters, indicating low metallicities in these systems. 
However, from our calculations we find that change of metallicities does not significantly af- 
fect the binary evolution when the CB disk is involved. The 11-min UCXB 4U 1820—30 is 
particularly interesting because of its negative period derivative P/P = —3.5 ± 1.5 x 10^^ 



^We notice two points about the persistent/transient criterion in UCXBs. Firstly, in Lasota et al. ( 2008h 
it is found that three of the five persistent UCXBs should be transient if their accretion disks are composed 
of pure helium or elements heavier than helium (C/0). However, in our CB disk model, the disks are not 
composed of pure helium but of mixed-compositions, so the accretion disks in these UCXBs are thermally 
stable, consistent with observations. Secondly, the transient source 4U1626— 67 should be persistent accord- 
ing to the criterion, but i t is transient in a different way: its outbursts do not last ten s of days but t ens o f 
years (jKrauss et al.l 120071 ). Th e value of Mi in T able 1 is calculated from Eqs. (4) in iKrauss et al.l (|2007l ) 
assuming a distance of 3 kpc (jChakrabartvlll998h and the time bet ween outbursts to b e 30 years. If the 
recurrence time is as long as 1000 yrs, it will yield a much lower Mi (jLasota et al.l 120081 ). and may help to 
resolve this transient/persistent problem. 



- 8 - 



yr~^ (Ivan der Klis et al.lll993al : IChou fc Grindlayll200ll ). which is inc onsistent with the lowe r 
hmit (-P/-P > 8.8 x 10~^ y^"^) of the standard evolutionary scenario (iRappaport et al.lll997l ). 
While previous explanations for this negative period derivative invoke acceleration of the bi- 
nary by a distant third companion in a hier a rchical triple system, or by the cluster pot ential 
han et all Il99ll : Ivan der Khs et all Il993bl : iKing et all 1 19931 : IChou fc Grindlavl boOlh . our 
CB disk scenario may present an alternative interpretation of this negative period derivative 
(see Fig. 5). 



4. Discussion and Conclusion 



During RLOF mass transfer, a CB disk may be formed as a result of mass outflow from 
the accretion disk, a nd has been invoked a s an efficient process for the removal of orbital 
angular momentum (ITaam fc SpruitI l200ll ). We propose a scenario for the formation of 
UCXBs from L/IMXBs with the aid of a CB disk in this work. The suitable binary parameter 
space (M2,i and Pi) with reasonable choice of the CB disk parameter 6 for the formation 
of UCXBs withi n 13.7 Gyr is found to be significantly larger than in previous "magnetic 
capture" model ( van der Sluys et al. 2005a| jbl). This difference is caused by the fact that the 
bifurcation period is considerably increased if the CB disk is included. In L/IMXB evolution 
the bifurcation period Pbif is defined as the initial orbital period when the donor star is 
on ZAMS, which separates the formation of converging systems from diverging systems. 
Because the value of Phif depends strong ly on the angular momentum loss mechanisms 
(Ivan der Sluys et al.ll2005al : iMa fc Li II2008I ). we would expect Pbif to be significantly changed 
when the CB disk is included in the binary model. In Table 2 we present the calculated 
values of Pbif for binary systems consisting of a 1.4 NS and a I.IMq secondary, with 
{6 = 0.002 — 0.01) and without a CB disk. We also list the corresponding v alues of the period 



Priof) at which RLOF begins (jPodsiadlowski et al 



2002 



Ma fc Li 1120081). According to the 



investigation of iPodsiadlowski et al.l (120021 ) and Ivan der Sluys et al.l (l2005a[ ). binary systems 
with initial orbital period slightly below the bifurcation period can achieve the shortest 
possible orbital period. We list the shortest periods for binary systems with certain value 
of S within 13.7 Gyr in Table 2, which clearly indicate that, (1) the larger 6, the larger P^f 
and Priof, and (2) UCXBs (P < 50 min) are not likely to form in such a scenario without 
the aid of CB disks (5 = 0). 

We note here that when the binaries reach their minimum periods where the donors 
become degenerate, their orbital periods will bounce back into the period-increasing phase. 
Our code cannot follow the binary evolution with a degenerate donor star. It is likely that 
the observed UCXBs may be explained as binaries both evolving to shorter orbital peri- 
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ods with a hydrogen-deficient, non-degenerate donor star (as presei ited in Figs. 6 and 7) 



and with a degenerate WD donor during a period-increasing phase (jYungelson et ahl 12002 



Nelson fc RappaportI l2003l : iDelove fc Bildsten] l2003l : iDelove et all |20oJ In addition, the 



hehum-don or channel (or semi- degenerate channel) may also contribute to the formation 
of UCXBs (ISavonije et al.l 119861 ) . Currently it is difficult to compare the (spectral) theo- 
retical models for hydrogen- deficient accretion disks with observations, or directly measure 
the orbital period derivative, one possible criterion to discriminate UCXBs in the period- 
decreasing/increasing phases is related to the donor masses, which, together with the orbital 
periods, allow to determine the mass-radius relation of the donor. Our calculations show 
that the UCXB's donor mass is around ~ O.IMq in the period-decreasing phase at P ~ 40 
mins, while in the period increasing phase, the d onor mass should be around O.OIM© (e.g. 
Yungelson et al.ll2002l: iNelson fc RappaportI l2003l ). Previously investigations on t he UCXBs 
XTE J1807 -294lFalanga et allboosh . SWIFT J1756.9-2508 JXrimm et al.ll2007h . and XTE 
J0929— 314 ( Galloway et al.l 20021 ) have shown that the donors should be WDs unless the 
binary orbital inclination is very small (< 10°), because of their small mass functions, while 
XTE 1751—305 is more likely to be in t he period-decreasing phase due to its relatively larger 
mass function (IMarkwardt et al.ll2002l ). Observationally there seem to be more systems in 
the period-increasing phase than in the period-decreasing phase. This may be a ddressed by 



the fact that UCXBs spend lon ger mass-transfer time in the former (> 10^ yr, seelRasio et al. 



2003 



2000; iDeloye &: BildstenI |2003| ) than in the latter phase (~ 10^ yr, see iNelson fc Rappaport 



and this wrok). From Table 1 there appears to be an apparent accumulation of systems 
with P ~ 40 — 50 mins. The reason may lie in that (1 ) during the period- incre asing phase, 
UCXBs spend more time at larger orbital periods (see lDeloye fc BildstenI |2003| . Fig. 9), and 
(2) when P > 50 mins, most of these systems bec ome transient sourc es with very weak accre- 
tion (< lO^^^M0yr~^). From the calculation of iDeloye et al.l (120051). these semi-de generate 



systems should mainly distribute at 40 mins < P < 90 min s (see lDeloye et al.ll2005l . Fig. 5). 



Most of them should be transients (see lDeloye et al.ll2005l . Fig. 2), among which the longer 



the orbital period, the weaker the accretion will be, thus the more difficult to be detected. 
These together account for the accumulation of UCXBs with P ~ 40 — 50 mins. Obviously a 
thorough population synthesis is needed to address the contribution to UCXBs from systems 
with non-degenerate, semi-degenerate and degenerate donors. 

However, t here exist sorne issue s in the CB disk scenario, especially the existence of CB 
disk in LMXBs. iDubus et al.l ( 12002| ) suggest that CB disks are prospective to be observed in 

'or th e exis- 
2008h 



infrared and sub-millimeter band. Althou gh observations have shown evidence 



tence of CB disks in young bi nary systems ( Monnier et al.ll2008t llreland fc Kraus 



netic cataclysmic variables (iHowell et al. 



200 



J 



Hoard et al.l 120071 ). and black hole LMXBs (iMuno fc Mauerhan 



Brinkworth et al. 



200 



i 



200 



i 



Dubus et a 



Gallo et al 



mag 



2007 



2003), 
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more infrared observations are still i ieeded to confirm or disprove the hypothesis that CB 
disk may exist in some LMXBs (e.g. iDubus et al.ll2004l ). Additionally, the CB disk param- 
eter 6 is poorly known, and it is possible to change with time or mass transfer rate. The 
strong dependence of LMXB evolution on the value of 6 prevents accurate estimation of the 
contribution of such binaries to UCXBs. More generally, we do not insist that there should 
be a CB disk in L/IMXBs, but argue that a mechanism that mimics its features may be an 
important ingredient for unde rstanding the overal l period distribution of UCXBs as well as 
cataclysmic variable binaries (jWillems et al.ll2005l ). 



Recent Chandra observations of nearby elliptical galaxies have revealed a population 
of luminous point X-ray sources , which are likely to be LMXBs with accretion rates M > 



10 yr ^ (e.g. lGilfanovll2004l : lKim &: Fabbiandl2004l ). These sources are explained either 



as transient LMXBs in which NSs accreting from a red-giant star in wide orbits (P > 10 d) 
( jPiro fc BildstenI I2OO21) or ultra-compact b inaries (P ~ 8 — 10 min) with a 0.06 — O.OSM© 
He or C/0 donor (IBildsten fc Deloyell2004j ). Our calculations suggest that normal LMXBs 
with a CB disk may present a plausible alternative interpretation for these luminous X- 
ray sources. We show in Fig. [8] the mean lifetime spent by LMXBs evolved to UCXBs at 
certain luminosity with S = 0.005. Here by the luminosity we mean the "potential maximum 
luminosity" , where the Eddington limit is removed and nearly all the mass lost by the donor 
is assumed to be accreted by the NS. From this figure we see that the UCXBs can be 
luminous {L > 10^^ ergs~^) for ~ 10^ yr, and the X-ray lifetime decreases sharply when 
L > 3 — 5 X 10^^ ergs~^, which may acc ount for the break in the luminosity function at 
~ 5 X 10"^^ ergs~^ (IKim fc Fabbiand 120041 ). A distinct feature of this explanation is that the 
luminous X-ray sources are predicted to be short-period, persistent rather transient sources. 
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Fig. 1. — Parameter space of the initial orbital period and donor star mass for binary 
systems which are able to evolve to UCXBs under the influence of a CB disk within 13.7 
Gyr. 
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Fig. 2. — The distribution of suitable CB disk parameter 5 and initial orbital period for 
progenitor binary systems with M2,i = I.IMq. 
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Fig. 3. — Evolutionary tracks of the donor star in a binary with M2,i = I.IM©, Pi — 1.04 day 
and different values of the CB disk parameter 5 in the H-R diagram. The circles and triangles 
indicate the beginning and end of mass transfer, respectively. The crosses correspond to the 
end of the calculation. 
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Fig. 4. — Evolution of the donor mass for a binary with M2,i = I.IMq, Pi = 1.04 day and 
different values of the CB disk parameter S. 
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Fig. 5. — Evolution of the orbital period for a binary with M2,i = I.IMq, P[ = 1.04 day 
and different values of the CB disk parameter S. The horizontal hue corresponds to P = 50 
mins. 
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Fig. 6. — Evolution of mass accretion rate (or half of the mass transfer rate) vs. orbital 
period for binary systems with M2 = I.IMq, Pi = 1.04 day and different values of the 
CB disk parameter 6. Th e red and blue line s indicate persistent and transient accretion 
according to the criteria in iLasota et al.l (120081 ) for a mixed- composition {X = 0.1, Y = 0.9) 
disk. Persistent (filled circles) and transient (triangles) LMXBs (including UCXBs) are also 
plotted for comparison. Here the open triangles mean that the derived mass accretion rates 
from observations are the upper limits. The symbols x, *, + on the evolutionary sequences 
denote where the composition of the donor in the binary is X = 0.3, 0.2, and 0.1 with 
Y = 0.98 -X and Z = 0.02, respectively. 
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Fig. 7. — Same as in Fig. 6 but for binary systems with M2,i = LlAf©, 5 = 5 x 10 ^, and 
p. = 0.94- 1.04 day. 
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Luminosity (erg s^) 

Fig. 8. — Mean X-ray lifetime spent by the LMXBs which could evolve to UCXBs within 
13.6 Gyr at certain luminosities with 5 — 0.005. The mean life time decreases sharply when 
L > 3 X 10^^ ergs-^ 
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Table 1: Parameters of known UCXBs 



Systems 


Type 


Porb (min) 


< M > (Mey-i) 


Refs 


4U 1820-30 


P 


11.4 


3.6tli X 10-9 


1 


4U 1543-624 


P 


18.2 


5.5+f X 10-10 


2 


4U 1850-087 


P 


20.6 


1.6 ±0.3 X 10-10 


3 


M15 X-2 


P 


22.6 


1.2 ±0.2 X 10-10 


4 


4U 1916-05 


P 


49.5 


8.9 ± 1.3 X 10-10 


5 


XTE J1807-294 


T 


40.07 


< 8±7 X 10-12 


6,10 


XTE J1751-305 


T 


42.42 


6±5 X 10-12 


6 


XTE J0929-314 


T 


43.6 


< 2± 1.5 X 10-11 


6,11 


4U 1626-67 


T 


41.4 


4.7+1^ X 10-10 


7,8 


SWIFT J1756.9-2508 


T 


54.7 


9.3 ±7 X 10-13 


9 



Note. — Here the type means persistent (P) or transient (T) UCXBs. The mean mass transfer rates 
< M > Usted here for persistent UCXBs are calculated by using the luminosities mentioned in the references 
and assuming accretion onto a IAMq NS with 10 km radius. The uncertainties in the calculated < M > 
mainly come from the uncertainties of source distance, neutron star ma sses and radii, spectral fit goodness, 
and recurrence times of transient sources (see iHeinke et al.l (|2007l |2009|) for a thorough discussion for these 
factors) . 
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Table 2: The bifurcation periods Pbif and Priofi and the shortest periods Pmin achieved for a 
binary with M2,i = I.IM© and different vahics of d. 



5 


^bif (day) 


^riof (day) 


Pmin (min) 





0.63 


0.5 


71 


0.002 


0.87 


0.75 


9 


0.004 


1.08 


0.92 


7 


0.006 


1.15 


0.96 


6 


0.008 


1.21 


1.04 


6 


0.010 


1.30 


1.13 


6 



